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ABSTRACT 

Many results in attitude analysis are still meaningful when the attitude is re- 
stricted to rotations about a single axis. Such a picture corresponds to attitude anal- 
ysis in the Euclidean plane. The present report formalizes the representation of atti- 
tude in the plane and applies it to some well-known problems. In particular we study 
the connection of the “additive” and “multiplicative” formulations of the differential 
corrector for the quaternion in its two-dimensional setting. 


Introduction 

I call our world Flalland, not because we call it so, but to make 
its nature clearer to you, . . . who are privileged to live in Space. 

— A. Square in Flatland 

The treatment of attitude, because of the non-linearity and non-commutivity of the compo- 
sition rule, is much more difficult to treat than position, for which components may be com- 
bined by simple addition. The complexity of the attitude composition rule leads to virtually all 
attitude problems being intrinsically three-dimensional or, in the case of the quaternion, four- 
dimensional. There is, however, a class of attitude problems which are much simpler, namely, 
single-axis problems, and the study of these will in many cases illuminate the more complex 
problems. The present report attempts to formalize such a treatment. 
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Attitude in Fiatland 


Having amused myself til a late hour with my favourite recre- 
ation of Geometry ; I had retired to rest with an unsolved prob- 
lem in my mind. 


Let us imagine that the world, Fiatland, has only two dimensions and a constant isotropic 
Euclidean metric. Such a world was imagined by Edwin Abbott Abbott [ 1 ], with the intent of 
satirizing the social and political foibles of his day as much as of clarifying the concepts related 
to the dimensionality of space. Our interest here is more limited than Abbott’s. We develop 
the mathematical structure of Fiatland somewhat further in order to better understand those 
aspects of attitude which do not depend on the dimensionality of space. The quotations which 
appear in this report are from [ 1 J. Following Abbott we will refer to our three-dimensional 
world as Space. 

In Fiatland, vectors are, of course, two-dimensional 


The “dot” product takes the usual form 

u • v = lijtq + u,v 2 , ( 2 ) 

while the “cross product” is now a scalar 

u x v = u 1 v. 2 - u 2 v [ . ( 3 ) 

There is, therefore, no vector product, and as alternate names to scalar and vector products we 
might prefer symmetric and asymmetric products. The lack of a meaningful vector product in 
two dimensions was ultimately the cause of many years of grief for Hamilton [ 2-4 J. 

The attitude matrix in two dimensions is a 2 x 2 proper orthogonal matrix, A, which trans- 
forms column vectors in the usual way 



W 

= A V , 

(4) 

with 






A r A = 

= aa t 

= 1 , 

(5) 


det 

A = + 

1 , 

(6) 

where 1 

1 denotes the 2 x 2 identity matrix, 





| 

'l 0" 




1 = 

0 1 


(7) 


It is a simple matter to show that in two dimensions the attitude matrix may be represented 
as 


cos # sin # 
-sin# cos# ’ 


( 8 ) 
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and 9 is the angle of rotation. If we define the matrix J according to 


J- 0 1 

J [-1 °J ’ 

(9) 

which satisfies 


j 2 = -i, 

(10) 

then Euler’s formula becomes simply 


A = cos 6 1 + sin 8 J , 

(11) 


which is much simpler than the three-dimensional form [ 5—10]. Note that J acting on a vector 
always generates a vector perpendicular to it. The matrices I and J in Flatland have an impor- 
tance similar to that of the 3 x 3 identity matrix and the Levi-Civita symbol in Space. They are, 
in fact, the representations of these objects in two dimensions. 

If we define now 

[[«]] = a J, (12) 


then trivially 

= (») 

which again is much simpler than the three-dimensional variant, and Euler’s formula becomes 

A — exp [[ 0 ]] , (I 4 ) 

as in Space. 

Corresponding to the quaternion in Space, in Flatland we must be content with the biemion 
(pronounced “by-Ernie-on” and named in honor ol Ernest P. Worrell, the character portrayed 
by Jim Varney). The biernion is defined as 


<1 = 


sin (0/2) 
cos(0/2) 


(15) 


for which 

q T q=l. (16) 

We continue to use the notation q (rather than b) in order to retain a greater resemblance to 
the equations in Space. 

In terms of the biernion Euler’s formula becomes 


A(q) = (<7? — ( i\) I + 2 q 2 <t\ J d 7a ) 

= {qj - q \ ) 1 + 2 q 2 [(<7i ]] (l ?b ) 

= (<7i I + ( h J) 2 • (l 7c ) 

The biernion may be extracted from the attitude matrix in a manner similar to the method for 

extracting the quaternion from the attitude matrix in Space, 

'/! = J V2 + tr/t , ( 1& >) 

(.4,, (18b) 
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where 


(19) 


trA — -f- A<>2 • 

Biernion composition follows directly from the trigonometric formula and reads 


q" = q’®q 

= { q' } 9 = { <1 } (}', 


(20a) 

(20b) 


where 




9 2 9 \ 

.-9l 9 2 - 


9 2 I + r/i J . 


(21) 


Note that biernion composition is commutative, as is the multiplication of attitude matrices in 
two dimensions. 

The Gibbs scalar or Rodrigues scalar is given by 


9 = 9i/ 9 2 = tan(0/2). 


Thus, 

9 

1J ’ 

and Cayley’s formula takes the familiar form 



, L+M 

i- W 

The composition of Gibbs scalars is given by 


in complete analogy to the formula for the Gibbs vector in Space. 
The Cayley -Klein parameters are 


a = ry 2 + i q x = e t6/2 , and 3 = q 2 - i q l = e“‘ e/2 = a c , 
and the superscript c denotes complex conjugation. These obviously satisfy 


It follows that 


a/3 = 1 . 

A = \ {or +/? 2 )l+ ^(a 2 -/3 2 )J. 


( 22 ) 

(23) 

(24) 

(25) 

(26) 

(27) 

(28) 


Attitude Kinematics in Flatland 

Restraining my impatience — for / was now under a strong 
temptation to rush blindly at my visitor and precipitate him into 
Space . . . 

The kinematic equation for the attitude matrix is given as usual by 

f t A {i) = [[u{t)]} A(t) , (29) 
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which, in fact, defines u>{t). If we define the biernion analogue, 




then the kinematic equation for the biernion is simply 

^-q{t)=\u{t)®q{i) = 7 ( 0 , 

at 2 2 

where 

= u> J . 

Likewise, we can partition { q } defined by equation (21) in terms of column matrices as 

{$} = [=($) g], 

which leads to 


^g(0= l=(g(0)w, 

at 2 


and 


H(g) 


<72 
- f /lJ 


Jr/. 


The kinematic equation for the Gibbs scalar becomes finally 

^ </(*) = 2 [1 + 9~(t)] u (t) ’ 


while that for the angle of rotation is just 


</0 

dt 


= u; . 


Euler’s equation for rigid-body dynamics is simply 


(30) 

( 31 ) 

(32) 

(33) 
(34a) 
(34b) 

(35) 

(36) 

(37) 


where N, the torque, is a scalar and I, the moment of inertia, another scalar, is given by 

/ = J r~ dm . (38) 


Attitude Errors in Flatland 

If Fog were non-existent, all lines would appear equally and 
indistinguishably clear. 

The representation of attitude errors in Flatland follows that in Space, with obvious simpli- 
fications. The error in the attitude matrix, since it has only a single degree of freedom, can be 
written as 

>1* = (MM true , (39) 
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with A * a random variable, usually an attitude estimate, and bA is the (multiplicative) attitude 
error, 

bA = exp{[[ A£ ]]} « I + [[ A£ ]] , (40) 

with A£, the attitude error angle, generally an infinitesimal quantity assumed to have zero 
mean. The attitude variance is defined to be 


P u = E{(At) 2 }, (41) 

where E{ ■ } denotes the expectation. 

The modeling of vector measurement errors follow a similar pattern. We write 


W = ell f I) /IV, 

(42) 

where t is a zero-mean random variable with variance a( v . In 
be written as 

linear approximation this may 

W = A V + AW , 

(43) 

with 


AW = f J AV = [[/lV]]e, 

(44) 

and we have defined [[ v ]] with vector argument to be 


([ v]] = Jv = [ V - . 

L- ", J 

(45) 


Thus, 


[[ u ]] 7 v = -u x v = -u 7 ’ [[v]], 

[[ v ]] T v = 0 , 

[[ u ]] T [[ V ]] = U V, 

[MHM] 7 ' = (u • v)l - vu 7 ', 

[[«]] V = [[v]]«. 

Batch Attitude Determination in Flatland 

I answer that though we cannot see angles, we can infer them, 
and this with great precision. 

We can now examine some well-known algorithms in their Flatland setting. These are the 
DYAD algorithm, the two-dimensional analogue of the TRIAD algorithm [11-12], and the 
BEST algorithm, the two-dimensional analogue of the QUEST algorithm [ 12]. The develop- 
ment of these algorithms in two dimensions is very similar to that of their forbears in Space. 
As can be expected, the results are much simpler in the smaller dimension. 


(46a) 

(46b) 

(46c) 

(46d) 

(46e) 
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The DYAD Algorithm 

For the DYAD algorithm we seek an attitude matrix which satisfies 

W = /IV, (47) 

where V and W are arbitrary vectors. In a space of n dimensions, n - 1 linearly independent 
vector measurements are required to uniquely determine the attitude matrix [13]. In two 
dimensions, therefore, a single measurement suffices. (In one dimension, zero measurements 
are sufficient.) 

To construct the attitude matrix we first construct orthonormal dyads of reference and ob- 
servation vectors as 


*1 

V 

= - 7 - and r., = Jr, , 
|V| 

(48a) 

and 



Sl 

w 

- _ _ and s., = Js, . 
|W| 

(48b) 

From 

J 3 = - J , 

(49) 

it follows that 

J A J T = A . 

(50) 

Hence, 

Sj = A r ; , i — 1,2, 

(51) 

Defining now orthogonal matrices 

(labeled by their columns) 


M h = [r 

i r 2 ], and M s = [s ( s., ] , 

(52) 

it follows that 

M s = A M n > 

(53) 

whence 

A = M s M'f l . 

(54) 


The development of the DYAD attitude variance follows almost identical steps as in the cal- 
culation of the TRIAD attitude covariance in Space | 12 ] with the result 


'Wad = <■ < 55 ) 


The BEST Algorithm 

The BEST (Bicrnion ESTimator) algorithm in Flatland is only slightly less complicated than 
the QUEST algorithm in Space. As usual, we seek an attitude matrix which minimizes [ 12, 

14 ] 

wn = ^ (56) 
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where the a t , i = 1, . . . , n > 2, are a set of positive weights, whose sum, we will assume, is 
unity. As in the Space we define a gain function, <j(A), such that 

g(A) = 1 - J{A) = tr (B r A ) , (57) 

which is a maximum when J(A) is a minimum, and, as before the attitude profile matrix is 
given by 

n 

B = (58) 

;=i 

Substituting equation (17a) in equation (57) leads straightforwardly to 

9(<i) = (<7> - q { ) * + 2 </, g, ^ , (59) 


where 


5 = tr (B 1 ) = tri? = B u + B 22 , 

^ = tr(J B t ) = -tr( J B) = B v2 - B 2l . 


Thus, 


</{<]) = <l' A'ry, 


with 



The maximization of <7(7) leads to the familiar eigenvalue problem 


(60a) 

(60b) 

(61) 

(62) 


K<l' = A inax <T, 

but in Flatland A niax can be calculated in closed form as 


(63) 




+ z* 


(64) 


and 


<1 = 


\R' + (^max + 


L \]»i + s . 


(65) 


The attitude variance of the BEST algorithm is calculated most easily from the Fisher informa- 
tion matrix using the fact that the BEST algorithm is a maximum-likelihood estimator [ 15]. 
Assuming the errors to have a Gaussian distribution, the calculation is straightforward and 
leads to 

U 

^iikst = y: <7 it , ~ • (66) 

t=i 


The optimal angle of rotation can also be computed directly by noting that the gain function 
can be written in the form 


(j{6) = s- cos 0 + z sin 0 , 


(67) 
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which is obviously a maximum when 6 = 6*, with 

cos e* = - * == and sin 6* = - = ■- . (68) 

V3 2 + Z 2 V5 J + Z 1 

We write the solution of equation (68) more conveniently as 

6 “ = arctan 2 (z,$) , (69) 

where arctan 2 is the same function as ATAN2 in FORTRAN. Equation (68) leads directly to 
a solution for the optimal attitude matrix, namely, 



Substitution of equation (65) into equation (17) leads somewhat less directly to the same re- 
sult, which should be compared with the construction of the optimal attitude matrix in Space 
developed by Markley [ 16], Markley’s FOAM algorithm ( 17] carries over with little change 
into Flatland and yields necessarily the same result as equation (70). 

General Comments on Attitude Estimation in Flatland ( and Space) 

I am about to appear very inconsistent. 

There seems to be some confusion concerning the use of representations in attitude estima- 
tion, which we will now attempt to muddy further. Typically in attitude determination, one is 
given a set of measurements, {zj , ... , z (V }, from which one wishes to infer the attitude, which 
we will denote without reference to a representation by A. The space of A we know from long 
experience is an m-dirnensional manifold, where m = 1 in Flatland, m = 3 in Space, and 
m — 6 in worlds so unfortunate as to be four-dimensional. An important milestone in ev- 
ery probabilistic estimate of the attitude is the construction of the probability density function 
(pdf) of the measurements as a function of the attitude, p xi |ZW ( z i» ... , z'^; A), where the 
primed variables denote the values taken by the (unprimed) random variables, and the attitude 
manifold is assumed to be simply a parameter space rather than a space of random variables. 
If A is also a random variable then the pdf of interest is p zi , ZN ,a( z ‘ t’ • • • » z 'n' A'). When 
one constructs a square loss function, one is, in tact, constructing part of the appropriate pdf 
assuming Gaussian random noise. 

The maximum-likelihood estimate is simply the value ot A (or A 1 ) which maximizes the 
appropriate pdf [ 18]. In mathematical notation we can write' 


^ML = al 'g max P Xl , 

(z!, 

... 1 ’ 

* z'n 1 A) 1 

(71a) 

A' ml = arg max p zi 

. , zpj , *4^1 ’ 

* 1 Z /V 5 ^ ) •> 

(71b) 


'In the particular case where A is a random variable one usually speaks of a maximum a posteriori estimate. 
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according to whether or not A is a random variable, and the maximum is taken over the man- 
ifold of A (or .4'). Likewise, in the special case that A is a random variable, we can define a 
minimum variance estimate of the attitude as 


A' m MV = E{A\z[,... ,z( v }, (72) 

that is, as the conditional expectation of A. This form of the minimum variance as given by 
equation (72) is not meaningful, generally, unless the representation of A has minimum di- 
mension. Otherwise, the conditional expectation will usually lead to a value which is not on 
the manifold, and, therefore, unacceptable as a solution. It is difficult, in general, to calculate 
the minimum variance estimate except in the case where the probability distribution of the 
measurement noise is Gaussian, in which case the minimum-variance estimate is identical to 
the maximum-likelihood estimate. 

The general method of solution by maximum-likelihood estimation to an attitude estimation 
problem given a set of measurements and a probabilistic measurement model is to write the 
negative-log-likelihood function 


J(A') = - log p %i ZNiA (z\, ... , z' N , A'), (73) 

where for definiteness we consider the case that A is a random variable. The negative-log- 
likelihood function is a minimum at the maximum-likelihood estimate. The procedure is 
thus to minimize the expression in equation (73) by an iterative method, such as the Newton- 
Raphson method. Thus, if .4' is the <-th iteration we write 


A' = AU3)0A'i, ( 74 ) 

where A{/3) denotes the general altitude as a function of /3, which is one of the many minimum- 
dimensional representations of the attitude which is Euclidean at the origin and for which .4(0) 
is the identity rotation. Expanding A as a function of / 3 leads to 


J(A') = J(A' i ) + 
1 


9(3 


J(A(P)®A' t ) 


0 = 0 


+ 2 ^ [d(3d(3 T ® A '>[ 


(3 + G( |/3| 3 ) , 


0 = 0 


(75) 


and minimizing this expression keeping terms only up to second order in (3 leads to the next 
iteration 


@i + 1 — 


' 0 2 
9(30(3 


jJiAWttA'i) 


-i 


0 — 0 


0 

dp 


J(AV3)®A'i) 


0=0 


A' l+ ,=A(P l+ } )®A',. 


(76) 

(77) 


This procedure will generally converge to a minimum of the ncgative-log-likelihood function. 
In well-defined attitude problems this minimum is usually unique and hence, 


A' m su = Jim A' t . 


(78) 
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In the limit that the amount of data is infinite the attitude covariance matrix can be written as 



d 2 

df3d{3 T 


J(A((3) ® A') 


l/3 = 0M' 


(79) 


If the measurement errors are Gaussian, then within the linearization approximation, equa- 
tion (79) will be true even for small samples. 


Treatment of the Biernion/Quaternion in Attitude Estimation 

It is high time that I should pass from these brief and discursive 
notes about things in Flatland to the central event . . . 

Several schemes have been proposed [ 19] for mechanizing the Kalman filter update for the 
quaternion. The effect of these and other schemes has been studied via numerical examples 
by Bar-Itzhack, Deutschman and Marklcy [20, 21 ]. These latter authors make a distinction 
between the update step of the Kalman filter using what they call the additive as opposed to 
the multiplicative update. This distinction is artificial and misleading, as we shall now show. 
Let us write the relation between the updated and predicted quaternions/biernions as 

?*(+) = «*(-) + Ar/,( + ), (80) 

which Bar-Itzhack et al. call the additive approach. The components are all resolved with re- 
spect to inertial axes. Let us examine the same equation expressed with respect to the predicted 
spacecraft body frame, i.e., we express all rotations as rotations from the predicted spacecraft 
body frame. Denoting the quaternions/biernions of rotation with respect to this frame by q' 
where 

<l'k = % ® ( 81 ) 

it follows that 

</!-(+) = 1 + Ar/J.( + ) , (82) 

where I = [ 0 0 0 1 ] T for quaternions, and 1 = [ 0 1 ] T for biernions. If we write now 

<*%-(+) = i ( 8 ^) 

then it follows that 

</*;( + ) = 7><7fc( + ) ® > ( 8 ^) 

which is the so-called multiplicative correction. Thus, the distinction between the additive and 
the multiplicative formulations of the Kalman filter is not one of the fundamental mechaniza- 
tion of the filter but simply the frame in which it is desired to compute the update. These 
two formulations are present in Reference [ 19], where they arc given the names “truncated 
covariance representation” and “body-fixed covariance representation.” Admittedly, the pre- 
sentation by those authors gave the appearance of there being one more distinct formulation 
of the Kalman filter than was actually the case. This has even led one careful study to test both 
formulations, as if they were distinct [ 22 ]. 

Where the important distinctions do lie is in how A q k or A q' k is calculated, and conse- 
quently, whether bq{-\-) has unit norm. From the earlier discussion it is clear that a correct 
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approach is obtained by expressing this quantity in terms of some representation of the atti- 
tude of minimal degree. In this case it is clearly advantageous to work from the spacecraft 
body frame so that this minimal-dimensional representation will be far from a singularity, and 
it will be most revealing to compare the results of References [ 19] and [20, 21 ] in that frame. 
The results of [ 20, 21 ], however, are not directly comparable to [ 19 ] because the former rests 
on the attitude Kalman filter of Bar-Ilzhack and Oshman [ 23 ]. However, many points of com- 
monality will be apparent. 

Consider now the estimation of a constant biernion from scalar measurements of the form 

z k = '$W k , k = l, . . . , N , (85) 

where u fc is a known direction in the spacecraft body and W k is some vector measured in the 
body frame. We assume that W k is related to a representation of the same vector in the inertial 
frame according to 

W k = AV k + v k , (86) 

where A is the attitude matrix and v k is white Gaussian noise. We wish to compute the batch 
attitude estimate from these measurements, using an good approximate estimate of the atti- 
tude as a point of departure. 

If we write now 

A = (SA) A a , (87) 

where A 0 is the approximate value of the optimal attitude estimate, then the measurement 
equation becomes 

z k = u r k (6A)W 0 ' k + uj:v k , (88) 

where W G k = A 0 V k> the expected value of the measurement in the body frame. 6 A is now an 
infinitesimal rotation, which we shall parameterize in terms of the additive biernion error as 
in equation (82). Recalling equation (17c) it is a simple manner to expand z k to lowest order 
in Aq with the result 

v k = Z k ~ Z o,k = 11 k A <7 + v k i ( 89 ) 

where k is the value of the measurement with A q - 0, v k is the scalar white Gaussian noise 
term appearing in equation (88) and the 1 x 2 sensitivity matrix H k is given by 

H k = [«i, k n 2 , k ] = [(u k xw 0<k ), (u,-W oJt )]. (90) 

The maximum likelihood estimate of A q (for the additive biernion correction, which is not 
constrained to preserve the norm) is given by 


^ qq P * 


(91) 


where the covariance matrix, P ' and the information vector, p, are given by 


p „ = 




Lfc=i 


-l 


N 


, p = E^:V 


(92) 


k-l 


For the multiplicative correction (which is norm-preserving) the estimate for the same data is 


, limit 1 nm!l V\ ^ 
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with 


-l 


^mult — 


N 


XX**r'"u 

* = 1 

V 


= (C)n ■ 




(93a) 


(93b) 


so that pi is just the first component of p. Note that we have written unnecessary (but not 
incorrect) transpose signs and not commuted symbols, even between scalars, to preserve the 
resemblance with the equations of Space. We can find a relation between the additive and the 
multiplicative corrections to the biernion by solving for p in terms of A</* dd and using the value 
of pj in equation (93). 

This leads to 

^9?, mult = A<7l,add + (P qq ) n {Pqq ) 12 A</2,add • $4) 

We will return to this equation soon. 

The additive correction, A^ dd , allows us to construct an optimal biernion, <?* dd . 


f /add - 1 + • (9^) 

Because it does not necessarily have unit norm, <y* dd does not without lurlher effort have an un- 
ambiguous connection to the attitude. However, we note that although </ add is not a biernion 
of rotation,” it is a sufficient statistic [ 18 ] for the attitude, certainly within the linear approx- 
imation of equation (89). It is, in fact, an estimate of the biernion of rotation, and we know 
also that were the measurement noise covariance to vanish (perfect measurements), g* dd would 
have unit norm and be the desired biernion. Thus, denoting the desired biernion of rotation 
by fi, we have that 

Cdd ~ r l + At7 at |d ’ (95) 

and 

A?) add ~ JV(0, /■*,,,) ■ (96) 

Hence, the negative log-likelihood function of </ add given // is 

d \v)=\ [(€, dd - n) T P qq (fold - n) + log del P qq + 4 log 2 tt] . (97) 

and the maximum-likelihood estimate of 7 / is simply 

T)' = arg max J I v), ( 98 ) 

7 ) : f) T ff~ 1 

where, since we know that the true biernion must lie on the manifold of unit four-vectors, we 
must maximize the negative log-likelihood subject to the norm constraint. 

We handle the constraint in the usual way, using Lagrange's method of multipliers, and 

optimize 

J(C Id I V) + 7> Xf l Tf l 
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without constraint and then choosing the value of the Lagrange multiplier, A, for which the 
constraint is satisfied. Differentiating the above expression with respect to 77 and setting the 
derivative equal to zero leads to 


*/* = (/ + APJ-’C,!, 


(99) 


and A is a solution of 


/(A) = i)* r (A) ?/(A) = <rj d (/ + XPJ- 2 £ dd = 1 . (100) 

We expect A P qq to be small. Therefore, it will usually be sufficient to calculate A using one 
iteration of the Newton Raphson method with vanishing initial value. Thus, 

^ ~ /'( 0) = 2 ('/add '/add ~ l) ■ (101) 

To first order 

'/add '/add - 2 -^'/_> • (102) 

Hence, 

* - ^'/. • ( 103 ) 

Substituting this in equation (99) leads to lowest order in Af/* dd 

r = + (104a) 

* Id (104b) 

= '/add “ ^'/ 2 ,add '/a’ 1,1 ■ ( 104c) 

The first component of the desired optimal biernion is simply (to this same order) 

7/1 = Ar/, add - {J } , IV ) 22 (/(,,, ) u Af/' a ,|d . (105) 

But 

- ( p Jn = (On (O12 7 (106) 

so that, in fact, comparing equation (106) with equation (94) wc have 

Vi = A(/[ (107) 

Since the other component must also agree to linear order in Ar/* lull , it follows that 

'I' = HnuU ■ (108) 

Thus, the additive correction to the biernion, followed by the normalization correction dictated 
unambiguously by the maximum likelihood criterion , is identical (at least up to linear terms in 
Ar/*) to the so-called multiplicative correction. It is hard to imagine that any other answer 
could be possible. It is obviously less burdensome to calculate the multiplicative correction 
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directly. Identical arguments hold for sequential correction of the biernion as in the Kalman 
filter. 

Discussion 

... my Lord has shewn me the intestines of all my countrymen 
in the Land of Two Dimensions . . . 

The representation of attitude in two dimensions has been described in detail. Two-dimen- 
sional analogues have been presented lor the well known TRIAD and QUEST algorithms. 
General properties of attitude estimation in two and three dimensions have been discussed. 
The question of whether the multiplicative or additive correction to the quaternion is to pre- 
ferred has been has a clear answer in Flatland. 

The additive correction, if done correctly, is identical to the multiplicative correction but is 
much more burdensome. The first commandment of biernion correction (and, one can show, 
also for quaternion correction in Space), thcrelore, is to multiply. We emphasize that this 
result is not the product of some heuristic argument or arbitrary procedure to be justified by 
experiment but the unavoidable conclusion to which one is led unambiguously and rigorously 
by the maximum likelihood criterion. 
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